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Abstract 

J\f =2 gauge theories broken down to = 1 by a tree level super- 
potential are necessarily at the points in the moduli space where the 
Seiberg-Witten curve factorizes. We find exact solution to the factor- 
ization problem of Seiberg-Witten curves associated with the breaking 
of the U{Nc) gauge group down to two factors U{Ni) x U{N2). The 
result is a function of three discrete parameters and two continuous 
ones. We find discrete identifications between various sets of parame- 
ters and comment on their relation to the global structure of = 1 
vacua and their various possible dual descriptions. In an appendix we 
show directly that integrality of periods leads to factorization. 



1 Introduction 

A breakthrough in our understanding of nonperturbative behaviour of A/" = 2 
supersymmetric gauge theories occured in 1994 with the work [1, 2] which 
showed that the low energy dyamics of the theory was encoded in the prop- 
erties of an associated hyperelliptic Seiberg-Witten curve. 

*e-niail: uf r j anikOif . uj . edu . pi 



1 



Recently there has been much renewed interest in this circle of problems 
associated with a set of completely new methods and insights into the struc- 
ture of A/" = 1 theories obtained by breaking M = 2 theories to A/" = 1 
by an inclusion of a superpotential for the adjoint A/" = 1 chiral superfield. 
Calabi-Yau constructions of the gauge theories, together with geometric tran- 
sitions [3, 4] gave rise to formulas for the exact superpotential which where 
later reinterpreted in terms of a computable random matrix theory [5] . This 
correspondence was subsequently proven in a purely field-theoretic context 
[6, 7]. 

Since the Af — 1 theories which arise in the above way sit in the region of 
Af — 2 moduli space where (some) monopoles become massless i.e. where the 
Seiberg-Witten curves factorize this new circle of ideas could be used as a new 
tool for studying the M = 2 phenomena and vice- versa. These interrelations 
were studied in [8, 9, 10, 11, 12, 13, 14]. In [15] matrix model methods were 
used to obtain a solution of the complete factorization of Seiberg-Witten 
curves for theories with fundamental flavours (see also [16] for an application 
of these results). 

The above dealt predominantly with the case when the gauge group was 
not broken. In the opposite case e.g. for the breaking U{Nc) — > U{Ni) x 
U{N2) it was found [17, 18] that the space ol M — 1 vacua exhibits a very 
complex structure of various connected components, each of which allows for 
multiple dual descriptions of the same physics but with different patterns of 
breaking. The analysis was later extended to the case of other gauge groups 
and matter fields in the fundamental representation [19, 20]. 

However in the preceding in order to map out the structure of the vacua 
one had to factorize the Seiberg-Witten curves for low Nc on a case by case 
basis. The aim of this paper is to provide an exact solution of the factorization 
problem for arbitrary A^^ and any pattern of breaking U{Nc) U{Ni) x 
U{N2). 

The plan of this paper is as follows. In section 2 we give a short intro- 
duction to the relation of factorized curves with J\f — 1 vacua. In section 
3 we summarize the equations for the meromorphic 1-form, the solution of 
which gives rise to the solution of the factorization problem. In section 4, 
which is the main part of the paper, we construct our solution. Then we 
proceed, in section 5, to study its properties relevant for the uncovering of 
the global structure oi Af — 1 vacua. We close the paper with a discussion 
and two appendices. In appendix A we give a mathematical self-contained 
proof that factorization follows from integrality of periods. We do this since 
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the condition of integrality of periods as a criterion for factorization was 
obtained using gauge-theoretical considerations. In appendix B we list the 
Mathematica code which implements our construction. 

2 Factorization of Seiberg-Witten curves and 
J\f = 1 vacua 

Since the celebrated work of [1, 2] it is known that N = 2 U{Nc) SYM 
theories develop an A^c-dimensional moduli space of vacua parametrized by 
vacuum expectation values of the adjoint scalar field 

(1) 

with 1 < p < Nc- The properties of the low energy dynamics of the theory 
around each such vacuum is encoded in the geometry of the Seiberg-Witten 
curve [1, 2, 21, 22] 

y2 = P^^(x,^,f -4A2^= (2) 

where the polynomial Pn^{x,Up) = (det(x/ — $)) = X^^=o " depends 
on the moduli UpS parameterizing the vacua through 

a 

OiSa + ^ kSa-kUk = (3) 
fe=0 

So = 1, Wo = (4) 

The H — 2 theory admits monopoles which at a generic point in the 
moduli space are massive. However at special submanifolds of the moduli 

space of vacua, n of these monopoles {n < Nc) may become massless and 
condense. At these points the Seiberg-Witten curve factorizes i.e. the rhs of 
(2) has n double roots and 2{Nc — n) single roots: 

PnA^, ul^'^f - 4A2^^ = • Hl{x). (5) 

The sets of moduli {m^"'^* } then depend generically on Nc — n continuous 
parameters and some discrete ones. 

The submanifolds where the monopoles condense are also important due 
to the fact that they are relevant when the Af — 2 theory is broken down to 
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J\f — 1 through an addition of a tree level superpotential 



n+1 



(6) 



Then the gauge symmetry U{Nc) may be either unbroken, broken down to 
U{Ni) X U{N2) with Ni + N2 = Nc etc. This is easy to understand clasically 
as it corresponds to redistributing Nc eigenvalues of $ among respectively 
one, two or more minima of the potential. On the quantum level the exact 
description is furnished by the factorized form of the Seiberg-Witten curve 
(5) with respectively n — Nc — 1, n — — 2 etc. 

The general solution for n = A^c ~ 1 (complete factorization) has been 
found in [23] . The solution is expressed in terms of Chebyshev polynomials 



where rj = exp^nik / Nc) . Note that this solution depends both on a continu- 
ous parameter xq and on a discrete one k — O...Nc — 1 (counting vacua in an 



Some solutions for low Nc in the case of breaking U {Nc) U (iVi) x U {N2) 
were found on a case by case basis in [17]. The aim of this paper is to provide 
a general solution in this case for any pattern of breaking characterized by 
A^i and N2 and the appropriate discrete parameters ki and ^2- 

Such a general solution is interesting especially as a very rich structure of 
the vacua was uncovered for small Nc in [17] (see also [18]). In particular the 
space of vacua turns out to have multiple connected components. Moreover 
within each component one may describe the same physics (the same Seiberg- 
Witten curve) in terms of different labels {Ni, N2, ki, ^2) and {N[, N2, k[, k'2) 
corresponding to different patterns of breaking of the gauge group. Note 
that one can define the A^j's and the appropriate glueball superfields S'j's 
not only in the semiclasssical limit (A 0) but also consistently at strong 
coupling through N^ = (l/27ri) £^ u where a; is a meromorphic 1-form (see 
next section) while ctj are some cycles. Dual descriptions correspond to 
different choices of the fundamental cycles a, but lead obviously to equivalent 
physics. The main motivation for this work was the fact that the knowledge 
of an exact solution of the factorization problem with different patterns of 

^The number of parameters is just Nc since the solution depends only on ry^. 





J\f^l U(Nc) theory^). 
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breaking of the gauge group might help to map out a global picture of the 
vacua and their possible dual descriptions. 

3 Integrality of periods 

A key object in the Seiberg-Witten theory is the meromorphic 1-form 

u;^T{x)dx^(tr^^ (8) 
which is exphcitly given for any Seiberg-Witten curve (2) through [8, 10, 7] 

00 = T(x)dx = , Pki^)'^^ ^ d /^^^^^ ^ _ ^^2^ X 

^P2^(x) -4A2^c dx \ V ) 

(9) 

From the definitions it is clear that 

1. T{x) has residue at infinity. 

2. The moduli of the vacuum Up can be reconstructed through 

11/ 1 
Up = -— ® -u: = -res2;=oo ■ ^ (10) 

3. The scale of the theory A can be reconstructed from the regularized 
integral 



oo 



reg 



uj\ = lim ( / u;- AT^logx 1 = -logA^^ (11) 



where a is a branch point of the Seiberg-Witten curve. 

4. If the Seiberg-Witten curve factorizes as in (5), the 1-form uj — T{x)dx 
defines a 1-form on the reduced curve 

y' = i^2(iv.-n)(a;) (12) 
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From the above properties we see that once we have found the meromor- 
phic 1-form for the factorized curve we can reconstruct the moduh Up and 
thus the solution to the factorization problem. We need therefore a criterion 
for the 1-form to come from a factorized curve. 

In [19] and in [18] (using different notation - see eqn. (A27) in the first 
paper of [18]) it was shown from the Konishi anomaly perspective and the 
Dijkgraaf-Vafa superpotential respectively that the periods of a; = T{x)dx 
should be integer^: 



The cycles a and (5 in (13)-(14) are ordinary compact cycles. As mentioned 
above uj can be understood to be a meromorphic 1-form on the elliptic curve 



The main aim of this paper is to give an explicit construction of the general 
solution to these equations. 

Note that the equations (13)-(14) were obtained from gauge-theoretical 
considerations. In appendix A we give a purely mathematical proof which 
shows that factorization indeed follows from the integrality of periods. That 
construction also shows a strong structural similarity between the standard 
Chebyshev polynomials appearing in (7) and the polynomials constructed 
here which solve the factorization problem. 

4 Construction of the meromorphic 1-form 

Before we proceed to solve equations (13)-(14) let us note some obvious sym- 
metries of the problem. Firstly, given a solution P{x), the shifted polynomial 
P{x) = P{x + xq) is also a solution. Secondly, a rescaling x ax corre- 
sponds to a rescaling A^^'^ —>■ A^^^/a^^". This allows us to change A and 
also generate distinct solutions for fixed A by picking a to be an appropriate 
root of unity. 

^We quote the equations for the case of n = Nc — 2. 




(13) 



(14) 



y"^ — F^i^x) = {x — a){x — h){x — c){x — d) 



(15) 
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Figure 1: The torus with modular parameter r and two marked points ai 
and 0,2 which will be mapped to infinity. 



The nontrivial part of the problem is then to find solutions parametrized 
by A^i and Ak and one continuous parameter. 

Let us note that in the original variables the problem seems intractable. 
Namely we have to solve the equations 



1 

TT J a 

1 r 

TT Jb 



NcX + II 



^/(x — a)(x — b) (x — c)(x — d) 



NcX + jj, 



\/ {x — a){x — h) {x — c){x — d) 



dx 



dx 



Ak 



(16) 
(17) 



for two of the five unknowns a, b, c, d and ii. These are coupled very nonlinear 
transcendental equations. And there seem to be no way to obtain an explicit 
solution from the above expressions. 



The torus with mcirked points 

The key to our solution is to use the representation of the elliptic curve (15) 
as a torus. However since the points at infinity play a special role in our 
setup we have to pick two points di and 5,2 on the torus (see fig. 1). These 
points will get mapped to the x = oo infinities on the two branches of the 
curve (15). 

The advantage of this representation is that holomorphic (and meromor- 
phic) 1-forms have a very simple form and periods may be calculated quite 
explicitly. Indeed using the general representation of a meromorphic function 
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with prescribed poles and residues on a torus [24] we may write lo as 



where 

oo 

,2 



Q[z) = e'''^'*'+2^^"^ (19) 

ri=— oo 

is the standard theta function and r is the modular parameter of the torus. 
Then cv has poles with residues Nc, —Nc at di — + ai and 0,2 — ^ + 02 
respectively. 

Now it is trivial to calculate the periods. Using the periodicity properties 
of the 9 function we get immediately 

— /'a; = — (20) 
27ri Jo 2ni ^ ' 

1 r Ct 

-y^. = iv.K -«.)+- (21) 

Of course due to the fact that lo has poles with residue ±A^c) the above 
expressions are only defined up to a multiple of N(.- We may now solve for 
the free parameters in a; to get 

C = 27riA^i (22) 

_ A^iT - Afe 
a2 — Oi = Aa = — 

'c 

In fact the final result will not depend on oi so we will set it to zero. 
Embedding of the torus 

In order to complete the solution of the factorization problem we need to 
represent the torus with the modular parameter r and the marked points 
ai = and = Aa as the curve (15). Once this is done we will have at 
our disposal the function x{z) which will allow us to reconstruct the moduli 
parameters Up and the scale of the theory A. Comparing the general forms 
of the meromorphic 1-form in both representations: 

^'"^ ^^ dx = {NJ{z) + C)dz (24) 

y 
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leads to the choice x — f{z) and y = dx/dz = f'{z). Although this choice 
is not unique, we will use the freedom of rescaling x — * ax + Xq at the end 
of the construction. In principle the knowledge of x{z) should be enough to 
calculate the Wp's through (10), however it is more convenient for numerical 
calculation to reconstruct explicitly the curve (15) and to calculate UpS from 

_ 1 p N^x + ii 

Up — res x=oo ^ ' , \^ziO J 

P ^y (x — a)(x — b)(x — c) (x — d) 

In order to reconstruct the equation satisfied hy y = f\z) and x = f{z) 
one can use standard reasoning and look at the Laurent expansion of f{z) 
around it's pole at z — di: 

f{z) K- + fo + fi{z- ai) + f2{z- d^f + h{z- (26) 

z — ai 

Then one fixes the coefficients of the quartic polynomial (15) by requiring 
that in the Laurent expansion of the difference 

f{zr - {f{z) + S,f{z) + S,nz) + S,f{z) + ^4) (27) 

all poles in [z — di) and the constant term cancel. Then the function (27) 
could have at most one pole at z — d2 but on a torus that is impossible. 
Therefore (27) is holomorphic hence constant (and equal to zero). The elliptic 
curve (15) corresponding to our solution is then 

y2 ^ + SiX^ + S2X^ + S3X + S4 (28) 

with the S'i's given by 

S, = -4/0 (29) 

-^2 = 6/0' -6/1 (30) 

Ss = 12/0/1-4/3-8/2 (31) 

-^4 = /o +7/1'- 6/0V1 + 8/0/2 -10/3 (32) 
The parameter /i in (25) is just /i — 2m Ni. 

The scale A of the theory 

It remains to evaluate the regularized integral related to A. It is convenient 
to calculate the equivalent quantity 



-logA^^^ = lim 
£^0 



a2—e 



uj + 2 log e 

ai+e 



(33) 
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The result is 



log = 21og^'(f) - log (^(f - Aa)9{f + Aa)) - — — - in (34) 

where f = (1 + t)/2. 

This completes our solution. We have constructed a solution parameter- 
ized by two discrete parameters Ni, Ak and the modular parameter of the 
elhptic curve r in the upper half plane. From this the standard rescaling 
X ^ ax + Xq provides the second continuous parameter xq, allows us to fix 
A to the desired physical value and introduces the third discrete parameter 
through the appropriate root of unity in a. 

Exphcitly we set 

a = -iA-^ ■ 9'{t) [Oir - Aa)9{f + Aa))"^ e"W • e"^"^ (35) 

where A is the physical scale of the theory, C — 2mNi, Aa = {Nit — Ak) /Nc 
and k is the new discrete parameter. 

In appendix B we give the Mathematica code which implements the above 
construction. 



5 Some properties of the solution and the 
global structure of Af = 1 vacua 

The factorization solution constructed in the previous section depends, for a 
given Nc, on the following set of parameters: 

[Ni,k,Ak,T,Xo] (36) 

where A^i, k and Ak are discrete, r lies in the upper half plane, while Xq is 
an arbitrary complex number. 

It may happen that two different sets of parameters describe the same fac- 
torized Seiberg-Witten curve. This may happen for two reasons. Firstly, as 
mentioned above the periods are defined only modulo so we may expect 
that there are transformation rules that allow us to reduce the parameters 
mod Nc- Secondly one expects e.g. that Ak has simple periodicity w.r.t. A^i 
and not only with N^. Thirdly, and this is the most interesting possibility, 
the different sets of (discrete) parameters for a single Seiberg-Witten curve 
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represent different dual descriptions of the same physics — this is part of the 
structure uncovered in [17]. Indeed the description of the identifications be- 
tween the sets of discrete parameters describes precisely the global structure 
of TV" = 1 vacua. 

If the Seiberg-Witten curve stays invariant, then necessarily the eUiptic 
curve also. Therefore identifications can only occur between sets of param- 
eters where r and r' are linked by a SL{2, Z) transformation. We will first 
consider identifications for fixed r and later deal with the cases t' = t + 1 
and t' = — 1/t. 



Identifications with t' — t 

The most trivial property is periodicity in k (see (35)): 

[TVi, A; + 27V„ A/c, r, Xq] = [N,, k, Ak, r, Xq] (37) 

Let us now consider a shift Ak — > Ak + Nc- Then Aa — > Aa — 1 and due 
to the invariance 6{z) = 6{z + 1) of the 6 functions the embdeding remains 
unchanged. The only modification is in a which gets rescaled by a phase 
^i-kNi/Nc^ Hence we have 

[TVi, k, Ak + AT,, T, xo] = [TVi, A; - TVi, AA;, r, xq] (38) 

The shift Ni ^ Ni + Nc is more complicated. Then C ^ C + 27riNc, 
Aa — > Aa -|- r. Using the transformation law 6{z + r) = e^'^'^e^'^^'''^6{z) we 
get X ^ X — 2'Ki and a ±a ■ p-^'^^^/^c^ The sign ambiguity is due to the 
branch cut of the square root in (35). Hence wc get 

[A^i + iVc, A;, AA;, r, xq] = [N^, k + AA;(+iVc), AA;, r, xq + 27ri] (39) 



Identification with r' = r + 1 

Under the modular transformation r ^ r + 1 the theta function transforms 
as 6{z, r + 1) = 6{z + 1/2, r). Therefore the elliptic curve remains the same, 
the only thing that changes is the identification of the periods namely 

^ N,{T+l)-Ak ^ N,T - {Ak - iVQ ^^^^ 

Hence we have 

[Ni, k, Ak, r + 1, xo] = [Ni, k, Ak - N^, r, Xo] (41) 
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This transformation law has a clear physical meaning. In our construction 
AA; = A-'i — k2 was periodic only modulo N^, while we know from the gauge 
theory perspective that ki, which labels vacua of U{Ni) should have peri- 
odicity A'^i. This identification is furnished by the modular transformation 
r ^ r + 1. 



Identification with t' — —1/t 

Here the modular transformation law of the theta function is more involved: 



0[w, 



(42) 



The meromorphic function defining x transforms as 



d , e(w) 
log 



dw e{w - Aa') 



T 



2m Aa' + — log 



e{z) 



dz 9{z - Aa'r) 



(43) 



where z — wt and Aa' = {—N[/r — Ak') /Nc- The periods may be calculated 
to be 

A^i = -AA;' AA; = N[ (44) 

Finally a transforms simply as a ^ ±t • a. This factor exactly cancels the r 
in (43) when constructing the torus embedding. Putting all these ingredeints 
together we get 



Ni,k,Ak, —,Xo 



-Afc,M+A^J,7V„r,^ + ^(^ + Afc 



(45) 



This is the key transformation which allows us to obtain a dual description 
with a different pattern of breaking U{Nc) U{Ni) x [/(A^a)- 



Global structure of vacua 

The above identifications allow us to study the global structure oi J\f — 1 
vacua. The connected components could be obtained by finding the orbits 
of the three discrete labels [Ni, k, Ak] under the group of transformations 
generated by (37), (38), (39), (41) and (45), assuming that there are no 
further identifications. We leave the detailed investigation of this structure 
to a subsequent work [25] . 
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6 Discussion 



In this paper we constructed an exact solution of the factorization problem of 
Seiberg-Witten curves with n = Nc — 2 massless monopoles. This is relevant 
for the breaking of jV = 2 theories down to M — 1 with the gauge group 
U (Nc) broken down to a product of two factors. 

The solution was obtained by constructing a meromorphic 1-form with 
integral periods on an elliptic curve. The sohition depends on three discrete 
parameters and two continuous ones. In principle one could generalize the 
structure to lower n by considering the construction of meromorphic functions 
on hyperelhptic curves of higher genus. Another direction in gencrahzation 
could be the extension to other gauge groups and/or addition of matter fields. 
A more complicated but interesting setting would be the consideration of 
quiver gauge theories, where nonhyperelliptic curves appear, and where it 
was shown that integrality of periods also plays an important role [26]. 

We also studied certain discrete identifications between the parameters 
and in appendix A we gave a proof that integrality of periods, a condition 
obtained in [19, 18] on physical grounds, indeed leads to a solution of the 
factorization problem. 

The physical interest of an exact solution lies in the possibihty of studying 
in detail the global structure of A/" = 1 vacua for any N(. along the lines of [17] . 
The main features of the analysis in [17], like the appearance of connected 
components of vacua and possible dual descriptions are related to discrete 
identifications between the parameters labeling the factorization solution. 
We leave a detailed investigation of this question for future work [25]. 

Acknowledgments RJ would like to thank the Niels Bohr Institute for 
hospitality while this work was carried out. This work was supported in part 
by KBN grants 2P03B09622 (2002-2004), 2P03B08225 (2003-2006) and by 
"MaPhySto" , the Center of Mathematical Physics and Stochastics financed 
by the National Danish Research Foundation. 

Appendix A. A mathematical proof of the fac- 
torization property of our solution 

The main object of this paper was to construct a solution of the factorization 
problem by constructing explicitly a meromorphic 1-form ou on an elliptic 
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curve 

,2 



y — {x — a){x — b){x — c){x — d) (46) 
with the following properties: 



J- (p cu eZ (47) 
OJ = ±^"e (48) 



These properties followed from the expressions for the superpotential in the 
Dijkgraaf-Vafa framework. Then we reconstructed the parameters of the 
factorized curve from the 'physical' formulas 

Up = -reSa;=oo x^ ■ uj (49) 
P 

In this appendix we show directly, without any physical input, that the 
conditions (47)- (48) give rise to a solution of the factorization problem, 
namely they allow us to construct a polynomial P{x) — x^" + . . . such that 
the Seiberg-Witten curve factorizes: 

p2 - 1 = (x - a){x - b){x - c){x - d) ■ H{xf (50) 

with H{x) a polynomial of degree Nc — 2. Prom the construction in the main 
text we see that we may assume without loss of generality that 

f u\ = \im ( f a; - A^clogxj = log2 (51) 

1 J reg \J a J 

Let US define the holomorphic function P{x) on the complex plane outside 
the cuts defined by (46) through 

P{x) = cosh (^J^ u?j (52) 

where a is a branching point of the curve (double covering) (46). Since the 
periods of u are integer this is a well-posed definition which does not depend 
on the path of integration. 

One can then show that P{x) (and its derivative) is continuous across 
the cuts. This follows from the fact that across the cut uj changes just by a 
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sign^. Hence P{x) defined by (52) can be extended to an entire function on 
the complex plane. 

Then the requirement that the residue of uj is Nc and the normalization 
integral (51) leads to the asymptotic behaviour 



J a 



cu Nclogx + \og2 + ... (53) 



which shows that P{x) must be a polynomial of order Nc with unit coefficient. 

We are now ready to show that factorization occurs. Let us denote F4 — 
{x — a){x — b) {x — c){x — d). Then the preceding footnote shows that we have 

p2 - 1 = F4 • G (54) 

We have to show that G has zeroes with even multiplicity. To this end 
consider the 1-form 

r- dP P'(x)dx 

but this is exactly equal to 

. — u sinh UJ I — , ^ 

V^- . ^ is = (56) 
smh / u> 

Ja 

which is a polynomial 1-form. Hence the roots of G have to appear in pairs 
which proves factorization 

p2 - 1 = F4 • (57) 

Note that the above reasoning works for any number of cuts and hence for 
any meromorphic 1-form with integer periods and fixed residues at infinity 
on an hyperelliptic curve. 

It is amusing to sec that the same construction generates the classical 
Chebyshev polynomials appearing in the solution of the complete factoriza- 
tion problem by Douglas and Shenker [23] . There, the relevant formula would 
be^ 

Nrdx 



PnA^) = cosh 



yx'^ — 1, 

The above formula indeed gives exactly the Chebyshev polynomial Tn^{x). 
Our formula (52) thus gives a very natural generalization of Chebyshev poly- 
nomials from the sphere to elliptic and hyperelliptic curves. 

''And consequently one can show that P is ±1 at the branch points of (46). 
"^Here for simpUcity we relaxed the normalization constraint (51). 
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Appendix B. Mathematica code which imple- 
ments the construction of the factorized Seiberg- 
Witten curve 

In this appendix we enclose the Mathematica code for generating the factor- 
ized Seiberg-Witten curve. The calhng sequence is 
generate [nc , nl , k , dk , tau , xO , dl] 

where dl stands for A^^=. All calculations have to be performed with high 
numerical accuracy. 

prec=30 (* numerical precision *) 

mer[a_,z_] :=Pi*EllipticThetaPrime [3,Pi(z-a) ,q] /EllipticTheta[3, Pi(z-a) ,q] 

(* torus embedding *) 

setpa:=({pa,pb,pc,pd,pe} = 

{(D[EllipticTheta[3, Pi(z-al) ,q] ,z] /. z->tt+al) , 
((l/2)D[EllipticTheta[3, Pi(z-al) ,q] ,z,z] /. z->tt+al) , 

((l/6)D[EllipticTheta[3, Pi(z-al) ,q] ,z,z,z] /. z->tt+al) , 
((l/24)D[EllipticTheta[3, Pi(z-al) ,q] ,z,z,z,z] /. z->tt+al) , 
((l/120)D[EllipticTheta[3, Pi(z-al) ,q] ,z,z,z,z,z] /. z->tt+al)}) 

setaa: = ({aa,bb,cc,dd} ={ pb/pa, -((pb''2 - 2*pa*pc)/pa''2) , 

(pb~3 - 3*pa*pb*pc + 3*pa~2*pd)/pa~3, 
-((pb~4 - 4*pa*pb"2*pc + 4*pa~2*pb*pd + 
2*pa"2*(pc"2 - 2*pa*pe))/pa"4) } ) 

seta:=({f0,fl,f2,f3} = 

{ (aa-mer [a2,z] ) /. z->tt+al, 
(bb-D[mer [a2,z] ,z] ) /. z->tt+al, 
(cc-(l/2)D[mer[a2,z] ,z,z] ) /. z->tt+al, 
(dd-(l/6)D[mer[a2,z] ,z,z,z] ) /. z->tt+al}) 

sets:= ({sl,s2,s3,s4} = {-4f0, 6(f0"2-fl), 12 fO fl -4f0~3-8f2, 
f0"4+7f l~2-6f0"2 fl+8f0 f2-10 f3}) 
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def curve :=setpa; setaa; seta; sets 
(* Seiberg-Witten curves *) 
u[i_] :=data[[i]] 
s[0]=l 

s [r_] : =N [-Sum [ s [r-a] a u[a] , {a,l,r} ]/r,prec] 
swpol [nc_] : =Suin [ s[a] x"(nc-a), {a,0,nc} ] 
swc [nc_,dl_] :=(swpol [nc] "2-4 dl) 
(* main procedure *) 

generate [nc_ , nl_ , k_ , dk_ , tau_ , xO_ , dl_] : =Block [{} , 
nu=nl/nc; deltak=dk/nc; 
al=0; 

a2=nu tau-deltak; 
q=Exp [Pi I tau] ; 
tt=(l+tau)/2; 

def curve ; 

Print ["Generating. . . "] ; 

thetaf ncs : =-I*Pi EllipticThetaPrime [3,Pi tt,q] * 
(EllipticTheta [3 , Pi (a2-al+tt ) , q] *EllipticTheta [3 , Pi (al-a2+tt ) , q] ) ^ (-1/2) * 
Exp[-2Pi I nl (nl tau-dk)/(2*nc''2)] ; 
al=dl~(-l/(2nc)) thetaf ncs Exp[-I*Pi*k/nc] ; 
f orminit=(x+2 Pi I nu) /Sqrt [x~4+sl x~3+ s2 x~2+s3 x +s4] ; 
f orm=N [al*f orminit /. x->al*x+xO, prec] ; 

data=Table [-nc Residue[form*x"i,{x,Inf inity}]/i, {i,l,nc}]; 
data=dat a/ (-Residue [form, {x, Infinity}] ) ; 
Chop [swc [nc , dl] ] ] 
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